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Abstract 

The integrable quantum group spl q (2, l)-invariant supersymmetric t—J model 
with open boundaries is studied via an analytic treatment of the Bethe equations. 
An su(2) feature is seen to hold for states at or close to half-filling. For these 
states the eigenvalues of the transfer matrix of the t—J model satisfy a set of su(2) 
functional relations. The finite-size corrections to the relevant eigenvalues, and 
thus the surface effect on the spin excitations, have been calculated analytically by 
solving the functional relations. 

1 Introduction 

The integrable supersymmetric t—J model of strongly correlated electrons has a long 
and interesting history (see, e.g., |[], 0, |], [5] and refs therein). More recently, the 
construction of the integrable version with open boundary conditions || |3| had to await 
the systematic development of boundary integrability f| |7|. For open boundaries, the 
integrable Hamiltonian reads 




L-l 

+ i sin 7 (m -n L )-\ sin 7 ^ (n>jS* +1 - SjU j+1 ) + H+ + H~, (1.1) 

i=i 

where the boundary fields Hf are dependent on two arbitrary parameters £± 0. The 
operators Cj± (cffi) are spin up or down annihilation (creation) operators. The Sj = 
(Sj, Sj, Sj) are spin operators and rij the occupation number of electrons at site j. The 
operator V = ITf=i(l ~ n j^ n jl) forbids the double occupancy of electrons at one lattice 
site. 
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We begin by recalling the essential ingredients underlying the integrability of the 
Hamiltonian QLip . The i?-matrix (see, e.g., |2|, HJ) is given by 
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(1.2) 



is the spectral parameter and 7 (or q = e l7 ) is the crossing parameter. This matrix is 
a trigonometric solution of the Yang-Baxter equation 



R 12 (u-v)R 13 (u)R 23 (v) = R 23 (v)R 13 (u)R 12 (u-v). 
On the other hand, the boundary X-matrices are given in terms of 



(1.3) 



sin£ 



( sin(£ + v) 

e iv sin(£ - v) 










(1.4) 



e sin(£ — v) j 



with [| 



K-(v) = K(v,£_), 

K + {v) = q- 1/2 K(-v + 7/2, «e+)M. 



(1.5) 
(1.6) 



The crossing matrix M is 



M 



( 1 \ 
q 2 




V 



(1.7) 



g 2 y 



The fT-matrices satisfy the boundary version of the Yang-Baxter equation 

i2 12 (u - v)K^{u)R 21 {u + u)-fr 2 -(u) = K^{v)R l2 {u + ^"(u^ifa - v). (1.8) 
Following Sklyanin 0, the open boundary condition transfer matrix is defined as |fj, §] 
2» = K+(v)U ac (v)K; d (v)U db l (-v), (1.9) 



with U ac (v) the monodromy matrix defined as the matrix product over the i?'s, 

< c \^)^SW<SW---<\ rfi W- (1.10) 
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Here the indices c and d in parentheses are in the quantum space C 3 x C 3 x ... x C 3 with 
indices a and b in the horizontal auxiliary space C 3 as usual. The operator £/ _1 (t>) is the 
inverse of U (v) , with 

U-J${v) = Btf 1 {v)R b t;f c 2 2 (v)Rlltl {v)...R a b d L L CL (v) , (1.11) 

where 

R*(v) = , -. Ud \ { V -. (1.12) 

cdy ' sm(7 + v ) sm(7 - v ) 

The elements of R will be denoted with a tilde, i.e., a ,b ,c± and w. 
Given the above, the commutation relations 

[T(v),T(u)) = (1.13) 

are fulfilled. The Hamiltonian H (|1 . 1|) follows from || [| 

dT(v) 



dv 



-- sinjH trK + (0) +trK + (0). (1.14) 



b=0 

In the following we consider the quantum group spl q (2, l)-invariant case only, for which 
Hf = 0. The transfer matrix ( |1.9|) has been diagonalised for this and the more general 
case by means of the algebraic Bethe Ansatz || [| . The eigenvalues are given by 



A(u) = Xa{v) + XvM + ^Vzx{v), (1.15) 



where 



^ / \ TT Cb(Vi — V )b(v i + V ) T / \ T / \ i / \ 

X Av) = UjT 1 w — l a L (v)a L (-v)k A (v), 1.16 

t=x K v i ~ + v) 

. / \ , , , 6(2^)6(2^ + 7) /c+(2v) \/ c„(2i; + 7)\ 
A Ct (v) = k VT (v)b L (v)b L (-v) , ( \ -A- / r - 1 U r ^ 

y 6(u - i;,-)^-!; - vt - 7) /ii 6(^i - u)a(i^ + v + 7) ' 
x / \ , / \ iL / \ 7l / v b(2v)b(2v + -f) fc+(2v) A 

x f 1 - + A n g0j^fa±I±7), (L18) 

V a{2v + 7) / 6(i/j - v)a(z/ j + v + 7) 

The related Bethe equations follow as 

/ a{v k )a{-v k ) \ L -it a(vi - v k )b(vi + v k )b(vj + v k + 7) 
V b(v k )b(-v k ) ) i-t k a{v k - Vi)a(i)i + v k )d(-Vi - v k - 7) 

'n ^rrlN ^ * = i jv, (i.i9) 

II IT 1 T = 1 > / = 1,-..,M. (1.20) 

i=i o(^ - Vi)b(-vi -Vi-j) 
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The local boundary factors are 

k A (v) = k Vl {v) = k Vxi {v) = 1. (1.21) 

The three possible states (|,|,0) represent either an electron (with spin up or down) or 
no electron (a hole). These are described by the numbers N and M of roots in the nested 
Bethe equations, where M is the number of holes and N — M is the number of down 
spins. As for the periodic case ||, it follows that the magnetization S z = |(nj — n\) = 
h(L — 2N + M) and the number of electrons Q = + n± = L — M are restricted to 
< S z < Q/2 < L/2. States in the half-filled band have one electron per site (M = 0). 

Our aim is to calculate the massless spin excitations of the states at or close to half- 
filling. This has been done for the periodic model in the isotropic (7 — > 0) limit via the 
root density approach to the Bethe equations |J. In Section 2 we establish an su(2) 
structure for the open quantum group invariant model and derive the corresponding 
functional relations. In Section 3 we derive the spin excitations by solving the functional 
relations for the finite-size corrections to the transfer matrix eigenspectra. The bulk and 
surface free energies of the vertex model and related t-J model are given in Section 4. A 
discussion and concluding remarks are given in Section 5. 



2 Functional relations 



We begin by considering the eigenvalue expression ( |1.15|) from the functional relation 
viewpoint. To show the su{2) structure we use semi-standard Young tableaux as in the 



study of the six-vertex model with open boundaries ]T0 



2.1 The su(2) structure 
Set C = 2L and define 

k 

|~T] = sin £ (7 — kj — v ) sin(2t> — 2j + 2kj) JJ 



N 



[~2~| = (?m(^ + ^7) sin £ (t> + A;7) sin(2t> + 2fc7)JJ 



sin(t> — v i + 7 + kj) sin(v + Vi + kj) 
\ sin(t> — Vi + k'j) sin(t> + — 7 + kj) 

sin(t> — Vi — 7 + Ary) sin(t>+t>j-|-/c7 — 2y) 
j sin(t>— Vi + kj) sin(t>+t>j — 7 + /C7) 



, (2.1) 



, (2.2) 



giafy+'Y+kj) si^ty—v—k'y) sin £ (f + /C7+7) sin (2v+2k / y—2 / y) sin(2f + 2/C7 + 27), (2.3) 



where 



^ sin{v - Vj + 7) sin(w + v 5 - 7) 
9m{v) = [I 



M sin(t> — Uj) sin(v + Uj — 27) 
The eigenvalue expression Q1.15 ) now reads 



(2.4) 



x/ n f^~? rn°\ sin- 1 (2i;-27) 
X(v) = (0 + ) — r7 - ■ \^ L/ L 



sin (7 + v) sin (7 — v) 



sin" 1 (2t; - 27) sin" 1 (2i; - 47) 
sin £ (27 — v) 



.(2.5) 
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To proceed further, introduce the auxiliary eigenvalues 



(i) 



T, 



E 



(2.6) 
(2.7) 



9+1 



where for each Young tableaux it is understood that there are relative shifts in the 
arguments: 



v+qy 








V 



(2.8) 

The zero superscript represents a shift in the right-most box. The number of terms in 
the sum is (q + 2), the dimension of the irreducible representations of su(2). Namely they 
are given by filling the numbers 1 and 2 in the (q + l)-box Young tableaux according to 
the rule that the numbers must not decrease moving to the right along the row. We thus 
get q + 2 numbered Young tableaux. We can show that the auxiliary eigenvalues satisfy 



or pictorially, 



rp(q)rp(l) _ jr(q+l) _j_ j- x J^ 1 ) 



□ 



TTT 



with 



(2.9) 



(2.10) 



fk 



The functional relations (|2.9f) can be further used to show that 



rp{q)rp{q) 

1 1 1 



k=0 



which coincides with ( [2.9|) for q = 1. It is also useful to introduce 

y « = T^TtV/llf k . 

k=0 



with Uq — 0. Then (|2.12j) can in turn be used to show that 



(2.11) 



(2.12) 



(2.13) 



(2.14) 



The relations (|2.9 ) and ( |2.12j ) are known as the T-system while (|2.14| ) is the |/-system 

ni pi p, m. 
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For half-filling (M = 0) we have S z = \L — N. Now gfjy(v) = 1 and the function 
is independent of M. It is quite clear from the su(2) functional relations (|2.9|) that the 
term /o contributes only to the bulk and surface free energies rather than to the finite- 
size corrections of higher order. However, the relation dT(0)/dv ~ dT^ (0) / dv implies 
that the finite-size corrections to the eigenvalues of the transfer matrix X(v) are governed 
by the auxiliary eigenvalues Tq. Thus the finite-size corrections to the Hamiltonian 
(1.1) follow from the consideration of T . As in the study of the six- vertex model with 



open boundaries 10], the surface effect on the finite-size corrections can be calculated 



analytically from the T-system (|2.9| ) and y-system ( 2.14 ) functional relations 



2.2 Zeros and poles 

The functional relations have been shown to be very useful in calculating the finite-size 



corrections to the transfer matrices of exactly solved models fl3| , 15, 1C ] . To solve the 



fusion hierarchy ( [2.9|) and (|2.14 ) we need to know the distribution of zeros and poles of 



the auxiliary eigenvalues and y( q \ We consider the model in the strip 

-7<Rew<7. (2.15) 

Inside this strip the transfer matrix T(v) in ( |1.9| ) is related to the super-symmetric 
t—J model ( |1 . 1| ) via ( p. .14 ). The clear advantage in working with the transfer matrix 



formulation is that it allows the application of powerful machinery from complex analysis. 
In this way we avoid the explicit manipulation of Bethe root densities, etc. The largest, 
or groundstate, eigenvalue of T^ 1 ' is not expected to possess zeros in the above strip. 
The zeros contributing to T are of order C from the bulk. Those contributed by the 
boundary are only of order 1, which become unimportant in the limit C — > oo. The bulk 
zero distribution is 

zero[T (1) (i;)] = 0, (2.16) 

q-2 

zero[T {q \v)] = \J{-k-f} c forg>l. (2.17) 

fc=0 

The zeros and poles of y"' are determined by ( |2.13| ), which gives 

(I) q = l : zero[t (1) (f)] = {0} £ , 

pole[tU(v)} = {- 7 } £ {7} £ , (2.18) 

(II) q>2 : zero[t iq \v)] = 0, 

pole^)] = {- g7 }%}£, (2.19) 

for the bulk contribution only. Here the boundary contribution to the zeros and poles, of 
order greatly less than L, are not listed and contribute less than those of order £ when 
the system size L becomes large. Only these zeros or poles of order L are especially 



important in the thermodynamic limit L — » oo [10 
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2.3 The functional relations for finite-size corrections 

The finite-size corrections to the eigenvalues can be obtained by solving the functional 
relations ( |2.9p and ( J2.14J ) in the strip ( 2.1 5[ ) . Denote the finite-size corrections of by 
TtL{v) and write 

T (1 \v) = TM e (v)Tlll(v). (2.20) 
The bulk and the surface free energy contributions together satisfy 

TB(v)T^(v + 1 ) = f . (2.21) 
Inserting ( ]2T20| ) into (p^2|) or (gj) we find that 

l£L(«)l£L(« + 7) = l + y (1) (v). (2.22) 

The finite-size corrections to (v) are thus represented by the ^-system component 
y^-'(v). In the following we give an analytical treatment of ( |2.22p and Q2.14| ). We will 
see that the finite-size corrections in the scaling limit are dependent only on the braid 
asymptotics and the bulk behavior of the functional relations. 

The Bethe equations ( |1.19|) -( |I.20|) render T^'(v) analytic. Since all functions involved 
in the eigenvalues are 7r-periodic, the analyticity domains for T^'(v) are not unique. It is 
thus useful to introduce functions of a real variable by restricting the eigenvalue functions 
to certain lines in the complex plane, 

r{x) ■= ^ ( l(^7 + ^), (2.23) 

a®(x) := y {q) (^xj + ^V), (2-24) 
A {q \x) := l + a {q \x). (2.25) 

For the groundstate the functions A^\x) and T(x) are analytic, non-zero (for those of 
order C) in —tt < Im x < tt and possess constant asymptotics for Re x — > ±00 (the 
ANZC property), which can be seen directly from the eigenvalues. 

Eqn ( |2.22|) can be solved using the new functions and applying Fourier transforms, 

1 roo 

F T (k) = — dx[\nT(x)]' e~ lkx , 

Z7T J-00 

[lnT(x)]' = / dkF T (k)e lkx } (2.26) 



in-, r i' 1 /.■• / 1 lh ' 

J —00 

A (q \k) = r dx [\nF ( /(x)}' e~ ikx 

271 J -00 



We obtain 



[\nF ( A q \x)}' = f°° dk A {q \k) e ikx . (2.27) 
J — 00 



lnT = A;*ln^ (1) , (2.28) 
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where the kernel k(x) is 



k(x):= r— — . (2.29) 

v ; 2vrcosh(2x) V ; 



Here the convolution / * g of two functions / and g is defined by 

/•oo roo 

(/ *g):= f(x - y)g(y) dy = g(x - y)f(y) dy . (2.30) 



There is an integration constant C in ( |2.28| ) which we drop because it does not contribute 
to the 1/L corrections. In case of the low- lying excitations we have to take care of zeros 
in the analyticity strips so that the simple ANZC properties hold. The result (|2.28|) is 
still correct if we change the integration path £ so that T(x) is ANZC and the Cauchy 



theorem can be applied as discussed elsewhere |L3, [10| . 



The function A^ is determined by the y-system (|2.14 ). According to Section [T2] 



the analyticity strip (|2.15 ) for y^ l \v) contains a zero of order C at u = and a pole of 



order C at u = ±7. All other functions y^ are analytic and non-zero in their analyticity 
strips —7 < u < 7. Taking care of these properties, applying Fourier transforms to the 
logarithmic derivative of the equations (|2.14 ) with the new functions ( |2.23| )-( |2.25| ) and 



then integrating the equations back we obtain the nonlinear integral equations 

lna (9) = lne (9) + k * InA^ + k * \nA ig+1) + D {q \ (2.31) 

where 

e (q) (x):=l 1 \ £n . (2.32) 
v ; \ tanh £ (|x) , q=l. v ; 

Here are integration constants. For the same reason we have to take care of the 
ANZC property in the analyticity strips in ( [2.31|) . 



2.4 The functional relations in the limit L —>■ 00 

The finite-size corrections can be extracted from the nonlinear integral equations ( |2.31| ) 
and ( p.28| ). The system size L enters the nonlinear equations (|2.31| ) through fl2.32Q . The 



function e^ 1 -* has three asymptotic regimes with transitions in scaling regimes when x is 
of the order of — In L or InL. We suppose that and A^ scale similarly. Thus in the 
following scaling limits, 



e { l\x) := lim e (q) ( ±(x + In C) 



L— >oo 



af (x) := lim a {q) ( ±(x + In £)) , (2.33) 
A { l\x) := lim A {q) (±(x + ]nC)) = 1 + a ( £(x), 



L— >oo 



eqn (|2.31|) takes the form 

h (g) = + k ^ + k + + D (q) ; (2.34) 
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where we use the abbreviations 



t {q \x) : 



x) := lna (9) (x) , £A {q) (x) := lnA (<?) (x), 

'0, q^l, 
-2e~ x , g^l, 



(2.35) 



and suppress the ± subscripts. The transfer matrix T(x) in the L — > oo limit now 
becomes 



lnT(x) = (A;*ln^ (1) )(x) 

1 r°° ( ln.4 (1) (?/ + ln£) ln„4 (1) (-y-ln£) 



27r 7-inA^lcosh(x — ?/ — ln£) ~^ cosh(a; + y + InC) ) ^ ~^ \C 

r e-m { *\y) dy+ e -^ T e^tA^\y) dy + </I 
Ln J-oo Lh j-oo \JL 

2 cosh x 
Ctt 



oo /I 

e-mW(y) dy + o -). 

-oo V JL, / 



(2.36) 



The above equation converges and can actually be evaluated explicitly with the help of 
the dilogarithmic function 



L(x) = - f dy — — — + ±ln:rm(l-:r). 
Jo y 



(2.37) 



Multiplying the derivative of ( ggg) with £A q , and ( ggf ) itself with (M 9 )', taking the 
difference, summing over g and using ( |2.35| ), we are able to obtain 



2jymt%)dy=-ZL(±_) 



(2.38) 



where the constants are given in terms of the x — > oo asymptotics by 



^2.39) 



The result ( |2.38|) shows that the finite-size corrections in the scaling limit depend only 
on the braid asymptotics and the bulk behaviour of the functional relations. 



2.5 Asymptotics and bulk behavior 

The nonlinear integral equations fl2.34|) can be easily solved in the limit x — > ±oo with 

7T 



7 = - h = 3,4, 



(2.40) 



In many cases different h correspond to different models. The equation Q2.38j ) shows that 
these asymptotic solutions are enough to obtain the finite-size corrections of the transfer 
matrix 
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It is obvious to see that the x — > oo asymptotics corresponds to the braid limit 
u — > ±ioo. In this limit (|2.14 ) reduces to 

(^ ) ) 2 = (l + ^ 1) )(l + ^ +1) )- (2-41) 

This equation in turn means 

2h iq) = £A {q - 1} + £A iq+1) + D {q \ (2.42) 

in terms of the functions a q . The constants D q can be either zero or non-zero as the 
different branches can be taken for the logarithmic functions in the nonlinear integral 
equations. 

To solve for y«) we write as 

v w = — (2 43) 

where the parameter 9 is to be determined. The recursion relation Q2.41| ) implies 

y { £ = sinq9 sin(g + 2)6>/sin 2 6», 
y® + l = sm 2 (q + 1)6 /sin 2 6, (2.44) 

for all q — 1, 2 • • -. This solution has to be consistent with the braid limit of T^(v). To 
fix the constant parameter 9 let us consider the groundstate, for which N = \{L + M) 
and 

lim (v) / <p(v) = 2cos(^J . (2.45) 

Inro-^±oo \hj 

Recalling Q2.43Q and using the relations 

y£ = lim yP/fo 

Imii- >±oo 



4cos 2 (tt//i) - 1 (2.46) 



hm r WTW// - 1 

lm!i->±oo 



we have 



= 7 = (2.47) 

Moreover, the special values of 9 lead to the closure condition 

yt 2) = ■ (2.48) 
For the sector S z = |(L + Af) — TV we have to modify to be 

777 77 

= m7 = (2.49) 

/ 6 
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where m = 2S Z + 1 = 1, 3, • • • < h - 1. 

In the limit x — > — oo, y( q > can be considered as the bulk behaviour in large L. Ac- 
cording to Section 1273 the analyticity strip for y^'(v) contains a zero of order N at v = 
and poles of order N at v = — ± 7. All other functions are analytic and non-zero 
in their analyticity strips in —7 < v < 7. For large A r we find that the leading bulk 
behaviour to 



constant 
constant 



taa(hhv) 



q = l 



(2.50) 



The constants are fixed by the functional equations (|2.14 ) and can be calculated similarly 
to the asymptotics of y( q \ As for the ABF model |T3, it is easy to see that the limit 



lim lim y 



(1) 



lim exp (— 2e 



. 



(2.51) 



Therefore the functional equations ( |2.14| ) are modified and we find the constants for 
2 < q < h — 3, with 



i/bulk 



sin(g — l)r sin(g + l)r/ sin 2 r, 
sin 2 gr/ sin 2 r, 



(2.52) 



where 



T 



m it 



(2.53) 



which is consistent with the closure condition ( |2.48| ). 

The transfer matrix eigenspectra has only one "quantum number" S z . There should 
be one free parameter between m and m . For the largest (groundstate) eigenvalue the 
appropriate choices are m = m = 1. The open boundary t—J system under consideration 
possesses spl q (2, 1) invariance. In the case of a fixed M the Bethe equations and transfer 
matrix eigenvalues are similar to those of the six- vertex model with open su q (2)- 
invariant boundaries. Thus we suppose that m = 1 as in the study of the XXZ-chain 



16, 17| and the related six-vertex model [10]. The low-lying excited states are given by 



m > 1. 

The solution y^ k ( v ) ls gi ven by 



Vhl\\A v )ybl\k{ v + 7) 



+ 

4 cos 2 T. 



(2.54) 



Thus we find lastly that 



±2 COS T 



tan(-/if) 



(2.55) 
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3 Finite-size corrections 



The finite-size corrections are only dependent on the braid and bulk limits. In these limits 
the functional relations are truncated and the summation in ( |2.38 ) can be replaced with 



-°° 9 =1 



Aw 



h-3 



\ DW£A W 



g=l 



(3.1) 



Recall that the constants DW are dependent on the branches of the dilogarithmic func- 
tions in the nonlinear integral equations. The appropriate choice yields the correct finite- 
size corrections. Simply taking =0 is consistent with the asymptotics solutions given 
in Section 2.5. To take nonzero 

DW 

we need to single out the appropriate branches from 
the asymptotic solutions of the equations, as has been shown for ABF models . 



A relevant dilogarithm identity has been established by Kirillov |T8 
functions 



V 



(9) (J, 



sin(g + 2)ip sm(q<p) 



sin 2 (<£>) 



with 



<P 



2 + r 



1 < b < n — 1, 1 < q < r, 



< j < r 



Consider the 



(3.2) 



(3.3) 



It is obvious that they represent the solutions of the asymptotic equations (|2.34|) with 
r = h — 2 or the solutions of ( |2.52 ) for the bulk behaviour with r = h — 3. Then we have 
the dilogarithm identity 



s(j, r) 



1 



9=1 

6 



+ yW(j,r) J 
3r 6j(j + 2) 



+ 6 j 



(3.4) 



2 + r 2 + r 

Now in terms of the dilogarithm function the finite-size corrections ( |2.36| ) are expressed 



as 



lnT(x) 



cosh x 
Cn 



h-3 
9=1 



Aw 



h-3 



+ \ D^lA^ 



9=1 



+4 



(3.5) 



Note that the nonlinear integral equations (|2.34|) including the closure condition ( |2.48| ) 



and their solutions presented in Section [T9are the same as those of the ABF models [13 



Therefore we can calculate the finite-size corrections in the same way. Similarly to 
with DW o, it can thus be shown that in terms of the functions s(j,r) the finite-size 
corrections (|3.5|) for the quantum-invariant open boundary system can be written 



In T(x) 



7r cosh x 
6C 



[s(0, h - 3) + s(0, 1) - 
-6(1 -m){2-m)\ 



s(m — 1, h — 2) 
1 

2 



(3.6) 



12 



Inserting ( |3.4j) into this equation we have the finite-size corrections in the more recogniz- 
able form [19| 

In T(x) = ^ (c - 24A m ) cosh x + e/~) , (3.7) 

where the central charge and conformal weights are given by 

6 



1 - 



h(h-l) 



_ [h-(h- i) m y _ i 

Am " Ah(h-l) ' ^ 



with m — 1, 3, • - • < h — 1. 



4 Free energies 

Consider now the bulk and surface free energies. Let 

t£!((v) 



Afreet) = T b (v)T s (v) 



sin(2w — 27) sin (7 + v) sin (7 — v) 



(4.9) 



where T b and T s are the bulk and surface contributions. They are determined by fl2.21| ). 
Factors of order L in fo contribute to T b and otherwise to T s . From (2J3) and (|2.21|) we 
should have 



T b {v)T b {v + 7) 
for the bulk and 

T s (v)T s (v + j) : 



sin L (7 + v) sin L (7 — v) 



sin 2L (7) 



sin 2i (7) sin L [7 + (7 + v)\ sin L [7 - (7 + v)\ 

sin(2w + 27) sin(27 - 2v) sin 2 (27) 



sin 2 (27) sin[7 + (7 — 2v)} sm[y — (7 — 2v)] 

for the surface. Solving these equations we find 

cosh(27A; — irk) sinh(w/c) sinh(7fc) cosh(t> k) 



\ogT b (v) 
\ogT s (v) 



2L 



dk- 



k cosh(7fc) sinh(7rfc) 



(4.10) 



(4.11) 



(4.12) 



00 cosh(47/c — irk) sinh(2t> k) sinh(27A; — 2vk) 
-00 k cosh(27/c) sinh(7r&;) 

cosh(27&; — nk) sinh(2t> k) smh^k — 2vk) 



dk 



k cosh(7/c) sinh(7rfc) 



(4.13) 



4.1 t—J model 



The eigenvalues E of the quantum invariant t-J Hamiltonian ( |1 . 1|) follow via the relation 
([LTD , with 



E~, 



4 d\{v] 



sin 7 dv [j 
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(4.14) 



From the above, we have that the groundstate energy is given by Eq 
where 



2L e b +2e s +ef+e-, 



2 d\ogT b (v] 



L sin 7 
4 



dv 



S1117 



v=0 

°° dk, C °^ 1 ^''^^ ~ sinh(7 k) 
-00 cosh(7fc) sinh(7rfc) 



(4.15) 



and 



2 dlogT s {v) 



+ 



S1117 
4 

sin 7 
4 



dv 



sin 7 



v=0 

°° dk C ° S ^^^ ~ sinh(27&;) 
-00 cosh(27&;) sinh(7r/c) 

°° <i£; COS k^^ _ sinh(7/c) 
-00 cosh(7fc) sinh(7rfc) 



(4.16) 



Here the local surface energies ef = for the quantum invariant case. At half-filling, the 
results for e& and e s are in agreement with those of the XXZ chain [16], as expected. The 
spin excitations are given by 



E 

Eq 



4 <91ogT flnit >; 



sin 7 



where the sound velocity is 



v s ir 



24L 



7T 



(c 



37 sin 7 

7T 



dv 
24A r 



7^0 
7 = 



v=0 

4' 



(4.17) 



(4.18) 



5 Conclusion and discussion 

We have exploited the su(2) structure of the transfer matrix functional relations to 
calculate the massless spin excitations of the integrable quantum group invariant t — J 
model (1.1) at close to half-filling. We took the special value 7 = 7r/h of the anisotropy 
parameter in order to close the functional relations. Results for the isotropic model are 
recovered in the limit h —>■ 00. Explicity, in this limit the scaling dimensions of the spin 
excitations, or spinons, follow from (|3.8|) as 



Xc 



2Ac 



z\2 



2(S Z ) 



(5.1) 



where S z 



0,1,... . 
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5.1 General parameter < 7 < 71/ 2 

The excitation spectrum can be obtained for general anisotropy parameter 7 via a related 
analytic nonlinear integral equation approach p0| , 
invariant model at or close to half-filling we expect 



log AO) = logT b + logT s 



6C\ 



c - 24A r 



2"B| , For the quantum 

1 



sm(— ) + 



(5.2) 



with the bulk and surface free energies given by (}4.12|) and ( [4.13|) . The central charge 
and conformal dimensions are 

6 7 2 



7r(7T — 7) ' 

[7m — (m — l)vr] 2 — 7 2 



(5.3) 

47r(7r- 7 ) ' (5 ' 4) 

where m = 2S Z + 1 = 1, 3, • • •. It follows that the finite-size corrections to the anisotropic 
t — J model are 



2Le b + 2e s - -^-{c-2AA 



(5.5) 



where the bulk and surface terms are given by Q4.15|) and (|4.16|) 



5.2 Local surface energies 

The above results are for vanishing boundary fields Hf. The quantum group spl q (2, 1) 
invariance is broken for Hf 7^ 0. However, the model remains integrable if H 

H~ = isin 7 (cot£_ - 1)^ - n h j2) (5.6) 
Ht = -isin 7 (cot£ + -l)(S£-n*/2) (5.7) 

where Hf = is recovered in the limit ^± — > 00. For finite £± the eigenvalues and Bethe 
equations are given by ( |1.15| )-( |IT20|) with || 

gsin(£ + — v) sin(£_ + v) 



sin £ + sin 

sin(£ + + v — 7) sin(£_ — v + 7) 
sin £ + sin 

sin(£_ — v + 7) sin(£ + + v — 7) 



(5.8) 
(5.9) 
(5.10) 



sin £ + sin 

At or close to half-filling the free energy can be calculated in a similar manner to that 



presented here for the quantum invariant case. As seen in Section 2.3, the total ground- 



state energy of the t-J model satisfies ( |2.21| ). The local surface free energy, on the other 
hand, satisfies 

q sin(£± + v) sin(£-|- — v) 



T±0)T±0 + 7 ) 



sin(f±) sin(f ± ) 



(5.11) 
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Here we omit the factor q in the following as it shifts logT s ± (f ) with no contribution to 
the surface energy of the quantum chain. We find 

j , f°° ^ efc(7_€±) cosh(2£ ± /c - irk) sinh(wfc) sinh(£±/c - v k) (512) 

j-00 k cosh(7fc) sinh(7rA;) 

Hence the local surface energy of the quantum chain is given by 

± 4 dlogT*^) 



trK + (0) sin 7 dv 



b=0 



4 r°° „ e fc(7 " 5±) cosh(2£ ± A; - vrA;) sinh(£±fc) 

are . (5.13) 



trK + (0) sin 7 J-00 cosh(7re) sinh(7rre) 

5.3 Hole excitations 

We have considered spin excitations, the so-called spinons, at or close to the half-filled 
band. The spinon excitations are given by the dimensions fl3.8| ) and ( [5.41 ), with ( |5.1[) 
for the isotropic model, where S z = \{L — 2N + M). One also needs to consider the 
holon part of the spectrum. According to the results obtained for the excitation spectra 
of the periodic t—J model, the holons should contribute another independent conformal 
theory ||. In contrast to the spinon part, the calculation of the free energies from 
X(v) is no longer as per Section |. Let us write X(v) = A froc (t> )A flnitcl (w)A flnitc2 (t> ), then 
Afinitei^) = Till e (v). Define A ircCiflnitc2 (w) = A free (t> )A finite2 (f ). Both the bulk and surface 
free energies should follow from 

A frC e,flnito 2 (^A froc , finitc2 (v + 7) = d(v)g M {v + 7), (5.14) 

which is clearly dependent on M, where 

= sin L (7 + y) sin L (7 - v) sin(27 + 2v) 
sin L (27 + 'y)sin L (-t;) sin(2v) 

This shows that the bulk and surface free energies follow from two parts, namely Qm(v) 
and d(u). The energies and e s given in ([4.15 ) and ( 4.16 ) are the contribution from d(v). 



The above relation also encodes the finite-size corrections to the eigenvalues, and thus the 
holon part of the spectra, through the contribution A flnitc2 (f). It is obvious that we need 
to analyse the Bethe equation (|1.19|) - (|1.20 ) to solve the inversion relation ( 5.14Q . It is 



possible that we will need to use a very different method to obtain the holon excitations. 
Whether or not the central charge associated with the holon part of the spectrum is also 
less than one for the anisotropic quantum invariant t—J model remains to be explored. 



Note Added. After completing this work we received a preprint by Asakawa and Suzuki 
in which the finite-size corrections are calculated for the open isotropic t—J model via the 
root density method |25| . They find the central charge c = 1 for both spinons and holons 
in agreement with the periodic case ||. The conformal weights of the spinons agree with 
our result Q5.1|) . In addition to treating the corresponding vertex model our results for 
the spinon conformal spectra generalise those of Asakawa and Suzuki to the anisotropic 
quantum invariant t—J model. 
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